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A limiting factor in the ability to interpret isotope efltect measurements in cuprates is the absence 
of sufficiently accurate data for the whole phase diagram; there is precise data for Tc, but not for 
the antiferromagnetic transition temperature Tjv. Extreme sensitivity of Tjv to small changes in the 
amount of oxygen in the sample is the major problem. This problem is solved here by using the novel 
compound (Cao.iLao.9)(Bai.65Lao.35)Cu30y for which there is a region where Tjv is independent of 
oxygen doping. Meticulous measurements of Tjv for samples with '^^O and find the absence 
of an oxygen isotope effect on Tjv with unprecedented accuracy. A possible interpretation of our 
finding and existing data is that isotope substitution affects the normal state charge carrier density. 



Isotope substitution is a powerful experimental tool 
used to investigate complex systems. Ideally, the iso- 
tope substitution affects only one parameter, for exam- 
ple, a phonon frequency which is directly related to the 
nuclear mass. However, the strong coupling of many pa- 
rameters in cuprates highly limit the ability to interpret 
isotope effect (IE) experiments. Eventhough the oxygen 
isotope substitution is known to affect the superconduc- 
tivity transition temperatures it is unclear whether it 
primarily impacts phonons, polarons, magnons, doping 
or other physical properties. :1] 

The isotope effect is usually described using the isotope 
coefficient a via the relation 
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where Tq is a phase transition temperature, M is the 
isotope mass and q = C, N, and g for the superconduct- 
ing (SC), antiferromagnetic (AFM) and spin glass critical 
temperatures, respectively. In many conventional super- 
conductors a was found to be very close to 0.5.0,01 The 
explanation of this a in terms of Cooper pairs glued by 
phonons was one of the triumphs of the BCS theory for 
metallic superconductors. [1| 

In cuprates the isotope effect is much more complicated 
and a is not single valued and varies across the phase 
diagram. The consensus today for YBCO like compounds 
is that in the SC phase, close to optimal doping, the 
oxygen IE is 
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which is very small but non zero.l5| On the SC dome ac 
increases as the doping decreases. In the glassy state 
less data is available, but it seems that the isotope ef- 
fect reverses sign and ag becomes negative. In extremely 
underdoped samples, where long range AFM order pre- 
vails at low temperatures, data is scarce, controversial, 
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FIG. 1. The (Cao.iLao.9)(Bai.65Lao.35)Cu30y phase diagram 
[l5l |. The antiferromagnetic, spin glass, and superconducting 
phases are represented in red, green and blue, respectively. 
For y < 6.6 the Neel temperature does not depend on y. This 
enables accurate measurements of the oxygen isotope effect 
on the Neel temperature Tjv. 



and has relatively large error bars.j^ The most re- 
cent measurements with Yj,Pri„j,Ba2Cu207_j show that 
aN = 0.02(3) in the parent compound. [5| There are also 
several theories dealing with the variation of a along the 
phase diagram, but since and aN are within an 

error bar of each other one cannot contrast these theo- 
ries with experiments. In particular, it is impossible to 
tell whether the same glue that holds the spins together 
holds the cooper pair together, or not. Increasing the ac- 
curacy of the IE measurements of the Neel temperature 
will shed light on the role of isotope substitution. 

As mentioned above, a major limitation on the mea- 
surements of aN is the strong dependence of Tn on 
doping. For example, in Yj,Pri_yBa2Cu307_5, Tn de- 
creases with increasing doping at a rate ATn = 2.5 K 
per Ay = 0.01. This strong temperature dependence is, 
of course, common to many other cuprates. As a conse- 
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quence, it is very difficult to prepare two samples with 
exactly the same T/v even with the same isotope. The 
smallest fluctuation in either y or d may lead to a huge 
fluctuation in Tjv regardless of the isotope effect. This 
is not the case for the (Cao.iLao.9)(Bai.65Lao.35)Cu30j, 
where T/v is constant for oxygen density y < 6.6, as can 
be seen in Fig. [Tifisj 

For our experiments, four sintered pellets were pre- 
pared using standard techniques, Two of the pellets 
were enriched with isotope and two with ^^O iso- 
tope in the same procedure: The samples were placed 
simultaneously in two closed tubes, each with different 
isotope gas, and then they were heated to allow the iso- 
tope to diffuse into the sample. In order to achieve a 
higher percentage of gas, the enrichment was repeated 
several times. 

The ^^O isotope content in the samples was determined 
based on measurements of gas composition being in equi- 
librium with the sample during the exchange. Balzers 
Prisma mass spectrometer was used to analyze in situ 
isotopic composition of the atmosphere. After the ex- 
change process was performed, the weight increase of the 
sample was also determined as the light ^^O isotope was 
exchanged with the heavy ^*0. The isotope enrichment 
in the samples measured by both methods looked to be 
higher than 80%. Finally a Thermal Analysis experi- 
ment (TA) was performed for the investigated samples 
after all experiments described in this work were fulfilled 
in order to verify the isotope enrichment. The samples 
were heated up to 1200°C in the NETSCH STA 449C 
Jupiter analyzer in a stream of helium. During the TA 
experiments, ion current signals for the ^^©2, ^^O^^O and 
^^©2 molecules were measured using a mass spectrometer 
(ThermoStar Pfeiffer Vacuum) . The results are shown on 
Fig.O The isotope content deduced from these measure- 
ments (comparing peak areas for the signals of particular 
oxygen molecules) was larger than 70%. 

The oxygen IE was measured using zero-field muon 
spin rotation/relaxation (/xSR). We particularly used the 
ISIS facility, which allows low muon relaxation rate and 
rotation frequency measurements. This is ideal for mea- 
surements near magnetic phase transitions where the 
muon signal varies on a long time scale. /iSR data at 
temperatures close to the phase transition are shown in 
Fig. 131 As the temperature is lowered from 383.4 K the 
relaxation rate increases. At T =378.5 K oscillations ap- 
pear in the data indicating the presence of long range 
magnetic order. The frequency of oscillations and the re- 
laxation rate increase as the temperature is further low- 
ered. The formula 

P^(t) ^ P^{ae-^'^ + {1 - a)e-^''*cos{ujt)) + Pne''^^ (3) 

was fitted to the muon polarization, where Pm, Ai , A2 and 
w are the polarization, relaxation rates, and frequency of 
muons spin in the fraction of the samples which is mag- 
netic, and a is the weighting factor between the muons 
experiencing transverse field and longitudinal field. This 
factor is close to 2/3 and temperature- independent. P„ 
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FIG. 2. Ion current signals for the ^^02, ^**0^'^0 and ^''02 
molecules (red, green and blue color, respectively) obtained 
during heating of the sample in a stream of helium. This 
graph shows that the isotope fraction in the samples is 
bigger than 70% (see text). 



and A are the polarization and relaxation of the spin of 
muons that stopped in the non-magnetic volume of the 
sample. The solid lines in Fig. [Sja) represent the fits. 

The AFM order parameter is determined by the fre- 
quency io = 7_B, where B is the average magnetic field 
at muon site and 7 is the muon gyromagnetic ratio. 
This frequency can be easily extracted from the /iSR 
data well below the transition but is very difficult to 
determined near the transition. A second approach is 
to treat the magnetic volume fraction Pm as the order 
parameter, [sl Il8l| P^ can be followed more closely to Tjv, 
although this parameter also has large error bars when 
uj is not well defined. We used an alternative approach 
similar to Ref. 19]; we define an order parameter which, 
does not require a fit, via the relation 



OV{T) 



^ < P,nf >-< P{T) > 

~ < P^nf >-< P{0) > 



(4) 



where < P(T) > is the average polarization at tempera- 
ture T, and < Pinf > is the average polarization above 
the transition. The denominator normalizes OV to 1 at 
zero temperature. All three order parameters for one 
sample are shown in Fig. [3fb) . The transition temper- 
atures determined using the different order parameters 
are in good agreement. The advantages of OV are clear: 
it is a model free and has very small uncertainties. 

In Fig. m we present the OV for two samples with ^^O 
and one with ^^O around 380 K. A wide temperature 
range from 50 to 410 K is shown in the inset for an ex- 
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FIG. 3. (a) The /xSR raw data at different temperatures and 
fits of Eq. [3] to the data in solid lines, (b) Comparison be- 
tween three different methods used to describe the AFM phase 
transition: (I) the muon precession frequency to obtained from 
Eq. [3l (II) the magnetic volume fraction Pm determined by 
Eq. O and (III) OP calculated by Eq. H 



periment done on separate occasion, in which two sam- 
ples of ^^O and one with ^^O were examined. We de- 
termine Tjv by fitting a straight line to the data in the 
main panel of Fig. 2] in the temperature range 378 to 
382 K, for each sample, and taking the crossing with the 
temperature axis. We find that T^8=382.49(0.34) K and 
r^6^382. 64(0.29) K. For 100% isotope substitution the 
isotope exponent is determined by 
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(5) 



When taking into account the isotopic fraction in the 
samples we obtain: 



ajv = 0.005 ±0.011. 



(6) 



This result indicates that < ct'c (see Eq. [5]) beyond 
the error bars, and is consistent with no isotope effect on 
Tn- 

One possible interpretation of these results is that mag- 
netic excitations are not relevant for superconductivity 
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FIG. 4. The parameter OP (see Eq. |4]) versus temperature 
near Tjv for and ^^O rich samples. The solid lines are fits 
to the data near the phase transision used to determine Tn. 
The inset shows a second experiment on the entire tempera- 
ture range. No oxygen isotope effect of Tn is observe within 
experimental error. 



since the isotopes affect Tc without affecting T/v. This 
approach was presented, for example, by Zhao et al..^ 
They found that samples enriched with ^^0 have longer 
penetration depth A than samples enriched with ^^0 with 
the same amount of oxygen per unit cell. A is related 
to the SC carrier density and effective mass m* by 
A~^ cx ns/m* so a priory both Us and m* can be affected 
by isotope substitution. ITf They ruled out the possibil- 
ity that the number of carriers concentration varies by 
demonstrating that the thermal expansion coefficient of 
samples with different isotopes are the same. The au- 
thors therefore concluded that the IE changes the mass 
of the cooper pairs, which could be explained by pola- 
ronic supercarriers. 

An alternative interpretation is that the isotopes affect 
the efficiency of doping, as suggested in Ref. [11] ■ To 
demonstrate this interpretation we present in Fig. [Sja) 
the critical temperatures in Y^Pri-yBagCuaOr-a for the 
two different isotopes taken from Ref. [5]. seems to 
be a bit higher than T]f, but T^^ is a bit lower than 
T^^. However, if we define an efficiency parameter Ki 
which relates the number of holes p to the number of 
oxygens in the unit cell y via p = Kiy, where i stands 
for the isotope type, we can generate a unified phase 
diagram. This is demonstrated in Fig. [Sljb) andlSjc) for 
T/v and Tc, respectively. In these graphs three different 
values of Kis are used to generate p while keeping Kie — 
1. When using Kis — 0.98, both curves of Tn and Tc 
versus p for the two different isotopes collapse to the same 
curve. Similar scaling of the doping axis was applied 
to the (Caa;Lai_a;)(Bai.75_a;Lao.25+a;)Cu30y system [ll| 
and was explained by NQR. [20| 
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FIG. 5. (a): Oxygen isotope effect of Tjv (left) and Tc (right) 
in YyPri_;,Ba2Cu307_i taken from Ref.E] The and ^"^O 
samples are in filled red circles and empty black squares, re- 
spectively. Figure (b) and (c) are demonstration of the IE 
using charge carriers density p = Kiy, insted of oxygen con- 
tent y. The values of Ki are given in the the figures. Green 
triangles {Ki = 0.98) represent reduction of 2% in the num- 
ber of charge carriers in the samples compared to the 
sample. In this case both Tn and Tc are functions of p 
regardless of the isotope. 



Next we discuss the IE on the stiffness in the above 
scenario. We assume that the effective mass of the SC 
charge carrier m* , the critical p where superconductivity 
starts Pcrit , and where Tc is optimal Popt , are not affected 
by the isotope substitution as suggested by Fig. [Sjc) and 
demonstrated in Ref. [lO] ■ The stiffness can be measured 
by the muon transverse relaxation rate cr, and it is ex- 
pected that: 17] 



C(P* - Pent) 



(7) 



where C is a constant. Dividing the differential of a from 
Eq. [7|by the relaxation at optimal doping yields: 
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dp 



Popt Pcrit 



Popt - P crit 



(8) 



The expected change in the stiffness due to isotope 
substitution can be calculated from Eq. |8] using popt — 1 , 
Pcrit = 0.42 (which are extracted from Fig. Efb)) and 
aopt = 3.0^is~\[21] Fory = 0.8, Ale = 1 and Xig = 0.98 
we get da = 0.083 pts^^ . This value is consistent with 
the experimental value of da — 0.08(1) /is~^ reported in 
Ref. [21|. In other words a 2% difference in the doping 
efficiency between the two isotopes can explain both the 
variations in phase diagram in Fig.[SlJa) and the variation 
in the stiffness. 

Our experiment shows that oxygen isotope substitu- 
tion does not affect the Neel temperature and therefore 
does not play a role in magnetic excitations. However, 
the isotope effect of Tc does not necessarily imply that 
phonons play a role in cuprate superconductivity. We 
show that an isotope dependent doping efficiency can ex- 
plain the variation in Tc and in the magnetic penetration 
depth between samples rich in ^^O or ^®0. We would 
like to thank the ISIS pulsed muon facility at Ruther- 
ford Appleton Laboratory, UK for excellent muon beam 
conditions. This work was funded by the Israeli Science 
Foundation and the Posnansky research fund in high tem- 
perature superconductivity. 
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